Identification of Extra Neutral Gauge Bosons at the LHC Using b- and t-Quarks 
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New Neutral Gauge Bosons, Z"s, are predicted by many models of physics beyond the Standard 
Electroweak Theory. It is possible that a Z' would be discovered early in the Large Hadron Collider 
program. The next step would be to measure its properties to identify the underlying theory that 
gave rise to the Z' . Heavy quarks have the unique property that they can be identified in the final 
states. In this letter we demonstrate that measuring Z decays to b- and t-quark final states can act 
as an effective means of discriminating between models with extra gauge bosons. 
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In the coming years, it is anticipated that the CERN 
Large Hadron Collider (LHC), a pp collider with centre 
of mass energy -y/s = 14 TeV, will reveal a new level 
of understanding of the fundamental interactions when 
it starts to explore the TeV energy regime. For a num- 
ber of reasons, including the quadratic sensitivity of the 
Higgs boson mass to radiative corrections, it is generally 
believed that the Standard Model (SM) is a low energy ef- 
fective limit of a more fundamental theory and numerous 
extensions of the SM have been proposed. Many of these 
extensions predict the existence of new neutral gauge 
bosons (Z ! ) and other s-channel resonances [l|, [2|, |3, |4J, [5| . 
If a kinematically accessible Z' exists, it is expected to 
be discovered very early in the LHC program. Once such 
an object is discovered, the immediate task would be to 
measure its properties and identify its origins. This is a 
difficult task and there is a vast literature on Z' observ- 
ables and analysis techniques. 

A key ingredient in determining the nature of a new 
resonance is to measure it's couplings to fermions. The Z' 
couplings to leptons can be measured using three observ- 
ables: the cross section to leptons, the forward backward 
asymmetry, Afb, and the width, Tz> [a]. For quarks, 
studies have shown that rapidity distributions can be 
used to separate u-quark couplings from <i-quark cou- 
plings [7|, [8j. However, these analyses are statistical in 
nature so there will always be contributions from the 
other type of quark. In contrast, the ability to identify b- 
and i-quarks in the final state can be a powerful tool to 
measure quark couplings that can be used to distinguish 
between models that give rise to Z' bosons. 

Previous studies have pointed out that third genera- 
tion fermions, top quarks in particular, can be used to 
search for extra gauge bosons [§, W\_ 11 L 13, LL3, LLJ, LL5| 
and to distinguish between models [13, [l6( ■ While some 
have noted the possibility of using third generation t- 



and 6-quarks to distin guis h between models of extra neu- 
tral gauge bosons [3, [l5|, [l7|, [l8|, this subject has not 
been fully explored. The ability to identify heavy quark 
flavours offers the unique opportunity to measure individ- 
ual quark couplings that is not possible for light quarks. 
In what follows, we describe a method of using b- and 
t-quark final states to distinguish between models of new 
physics that predict extra neutral gauge bosons 
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The primary challenges in these measurements will be 
the identification efficiencies for top and bottom quarks 
needed to make statistically meaningful measurements 
and the discrimination of the t's and b's coming from Z' 
decays from SM QCD backgrounds. 

To distinguish between models, we propose to use the 
cross sections a(pp — > Z' — ► bb) and a(pp — ► Z' — ► U), as 
described by the Drell-Yan cross section with the addi- 
tion of a Z' [fj, [2fJ at the LHC. We computed the cross 
sections using Monte-Carlo phase space integration with 
weighted events, imposing a rapidity cut on the final 
state particles of \r/\ < 2.5 to take into account detec- 
tor acceptances. We also included pr and invariant mass 
distribution cuts with values chosen to reduce QCD back- 
grounds as described below. In our numerical results we 
take a = 1/128.9, sin 2 9 W = 0.231, M z = 91.188 GeV, 
T z = 2.495 GeV and m t = 172.5 GeV [211 We use 
the CTEQ6M parton distribution functions [22j and in- 
cluded a K-factor to account for NLO QCD corrections 
[23] while NNLO are not numerically important to our 
results [3,j2F 
important [26| 



|. Final state QED radiation effects are 
but require a detailed detector level sim- 
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ulation that is beyond the scope of the present analysis. 
The Z' widths only include decays to standard model 
fermions. NLO QCD and electroweak radiative correc- 
tions were included in the width calculations [271 ]. 

An important challenge for this analysis will be to 
achieve sufficiently high b- and t-quark identification ef- 
ficiencies to provide the statistics needed to distinguish 
between models. The ATLAS and CMS collaborations 
have worked hard at estimating these values but expe- 
rience with real data will be required to obtain reliable 
values. We therefore present results for two sets of values, 



distilled from the literature, that we expect to bound the 
values that will eventually be achieved by the LHC col- 
laborations. Once the LHC experiments start to collect 
data, these values should be refined as experimeters gain 
experience and a better understanding of their detectors. 

For b identification efficiency, the ATLAS TDR gives a 
value of eb = 60% for low luminosity running and 50% for 
high luminosity running with 100 to 1 rejection against 
light and c-jets |28(. We will use the latter value which 
is appropriate to the high luminosities we assume. The 
rejection of fakes arising from light and c-jets can be im- 
proved considerably by requiring that both the b and b 
are seen. We therefore consider two cases for tagging bb 
events; 50% when only one b is observed and e b i = 25% 
when both the b and b are detected, independent of the 
dijet mass. Note that the bb detection efficiency is likely 
to be higher than simply using e^. 

The understanding of i-quark identification efficiencies 
is evolving. The top quark almost always decays into a 
6-quark and a W + boson (t — > W + b) with the W's sub- 
sequently decaying either into two leptons (eu e , fiu^ or 
tv t ) or into a light quark- antiquark pair (ud, cs) that in 
turn hadronizes. The single lepton plus jets final state, 
where one W decays leptonically and the other W decays 
hadronically, tt -> WWbb -> (lv)(jj)bb, has a BR- 30% 
of all tt events and is generally viewed as giving the 
best signal-to-background ratio. With suitable kinematic 
cuts and including the BR to (lv)(jj)bb, a recent ATLAS 
study estimates e t t ~ 4% [14|. However, reconstructing 
the invariant mass of the tt system will reduce this num- 
ber Q. The ATLAS TDR is slightly more optimistic, 
claiming the efficiency for detecting a M t j = 2 TeV res- 
onance of about 5% including the semileptonic mode 
BR while a CMS simulation obtains the lower value of 
e ti ~ 2% H|. Baur and Orr [H, |jj(J found that the t- 
quark identification efficiencies for this channel can be 
improved by using 2-jet and 3-jet final states with b-tags. 
The fully hadronic modes have a combined BR — 45% 
so utilizing the hadronic modes has the potential of im- 
proving the tt identification efficiency significantly. A 
method has been suggested to distinguish top jets from 
standard model backgrounds u sing substructure of the 



top jet 3l|, |32j. Kaplan et al. [3l| estimated that high 
Pt dijets can be rejected with an efficiency of — 99.99% 
while retaining — 10% of the tt pairs. By combining the 
different top decay channels and identification strategies 
it should be possible to increase the overall tt identifica- 
tion efficiency. Given that the subject of t-quark identifi- 
cation at the LHC continues to evolve, we assume a wide 
range of values of e t i, taking 1% and 10% for the low and 
high efficiency scenarios respectively. 

Another challenge for making these measurements will 
be to distinguish the Z' signal from the large SM QCD 
backgrounds. The invariant mass distribution for bb final 
states is shown in Fig. Q] for the SM QCD background 
and the signal for a Z' with a mass of 2 TeV for sev- 
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FIG. 1: Invariant mass distributions for the Drell-Yan process 
pp — > bb including a Z' with mass M z i = 2 TeV and the bb 
QCD backgrounds. The sets of curves correspond to E@(i{>) 
[lj, Left Right symmetric (qr/ql = 1) [331] , Simplest Little 
Higgs 1341. [3 511 . 3-3-1 model [3y|, and TC models (tan# = 
0.577) 9, 37]. A kinematic cut of Pt > 50 GeV was imposed 
on the 6-quarks. 



eral representative models. The QCD backgrounds were 
calculated using the WHiZard package [38| with O'Mega 
matrix elenrcnt generation [39] and as an independent 
check we also calculated the QCD cross sections using a 
simple Monte Carlo event generator with tree level ma- 
trix elements. We use LO QCD cross sections in our 
background calculations. While it is known that higher 

mmi, nlo 



order QCD corrections can be substantial 
corrections are highly dependent on the region of phase 
space being studied. As a crude estimate of the impor- 
tance of NLO correction on our results, we rescaled the 
LO QCD backgrounds by a factor of 1.4 and found this 
to have little impact on our results. 

The px distributions are quite different for the sig- 
nal and backgrounds with quarks coming from Z' de- 
cays having a much harder distribution than the back- 
ground events. The background can be reduced consid- 
erably by imposing a transverse momentum cut on the 
reconstructed final state t and 6's at some expense to 
the signal. The pt cut was varied and it was found 
that the optimum cut is approximately px Q > 0.3Mz' 
which reduces the background significantly compared to 
the signal. A stronger cut improves the signal to back- 
ground ratio but decreases the total signal and therefore 
increases the statistical uncertainty. The invariant mass 
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certainties originating in the parton distribution func- 
tions. The ratios are defined by 
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FIG. 2: Invariant mass distributions for the Drell-Yan process 
pp — * bb including a Z' with mass M Z ' — 2 TeV and the bb 
QCD backgrounds including a kinematic cut of Pt > 0.3M z i 
on the b-quarks. 



distribution for the signal and background are shown in 
Fig. Rafter applying the cut. 

The QCD backgrounds can be further reduced by con- 
straining the invariant mass of the final state fermions 
to \Mft — Mz>\ < 2.5Tz>- The window was chosen to 
balance the total signal against the signal to background 
ratio. We examined the model independent choice of 
\Mfj-Mz>\ < 0.07 M z >, but found that our results were 
not very sensitive to the precise choice of M i j window. 

Fakes from gluon-, light quark-, and c-jets are poten- 
tially problematic but there is a tradeoff between heavy 
quark identification efficiencies and mistagging that re- 
quires detailed detector simulations. Likewise, we defer 
detector resolution effects to more detailed future stud- 
ies. Other non-QCD SM backgrounds include Wbb+jets, 
(Wb+Wb), W+jets, etc final states. Baur and Orr have 
shown that these can be controlled by constraining the 
cluster transverse mass and invariant mass of outgoing 
jets (and leptons) to be close to m t [13| |3CJ • 

In addition to the QCD backgrounds and the ques- 
tion of heavy quark identification efficiencies, there are 
additional theoretical uncertainties in the cross sections; 
higher order QCD and EW corrections to the cross sec- 
tions, both initial and final state contributions, and un- 
certainties in the parton distribution functions. We can 
reduce some of these uncertainties by using ratios of 
heavy quark production to fi + ^ production; Rb/^ and 
Rt/fi- In particular, these ratios nearly eliminate the un- 
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where K q is a constant depending on the QCD and EW 
correction factors, and the factor of 3 is due to summation 
over color final states. Each of these ratios depends on 
only four couplings from each model. An analysis based 
on the location of a measured Z' in the R},/^ ~ Rt/ ^ pa- 
rameter space provides a means of distinguishing between 
models. 

We assume that a Z' has been discovered and it's mass 
and width measured [6|, [20( so that the appropriate Mqq 
cuts described above can be applied. It is expected that 
a Z' with Mz> < 2 TeV can be discovered early in the 
LHC program with approximately 1-10 fb _1 of integrated 
luminosity depending on the specific model. 

To obtain our results we calculate the expected num- 
ber of events and statistical error for signal plus back- 
ground for a given integrated luminosity and particle 
identification efficiencies, e„+„- , e^i, and e t j. The ex- 
pected number of SM QCD and electroweak events were 
similarly calculated and subtracted from the signal plus 
background events to give the predicted number of sig- 
nal events. From these intermediate results we obtained 
the ratios given in equations [T] and [2] with the errors cal- 
culated in the usual way by including both signal and 
background. We did not include uncertainties coming 
from luminosity and identification efficiencies. In the lat- 
ter case there is simply too big a range to include in an 
error, rather we show results for the two cases discussed 
above. 

Our results for R b /„ and R t /„ are shown in Fig. [3] for 
Mz' = 2 TeV. Fig.[3]Ja) shows results for the high fermion 
identification efficiency values with la statistical errors 
based on an integrated luminosity of L = 100 fb _1 . The 
low 6tf case would require higher integrated luminosity 
to distinguish between models so in Fig. EIb) we show 
statistical errors based on L = 300 fb _1 . The errors scale 
as 1/yX and very roughly like 1/ \/£qq so one can esti- 
mate how the errors will change with different integrated 
luminosities and heavy quark identification efficiencies. 

It is clear that most models can be differentiated using 
heavy quark final states. However some models such as 
the Eq(^) and SU(3) x U(l) anomaly free Little Higgs 
model give similar ratios so one would need additional 
input such as leptonic observables to distinguish between 
them. 

In summary, we demonstrated that, in principle, the 
decay of a Z' boson into third generation quarks can be 
used to distinguish between models of physics beyond the 
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FIG. 3: R b/li versus R t/fJ , for M z > = 2 TeV for the E 6 (x), 
E 6 (tp)^ 6 (n) [j., Left-Right Symmetric Model (g R /g L = 1) 
(LR) (H, Alternate Left-Right Model (rrJol = 1) (ALR) 
(HI, Simplist Little Higgs Model (SLH) [MEl, Littlest Higgs 
Model (cot6 H = 1) (LH) [i|ppL3-3-l 2U1D Model [H, TC 
- Topcolour (tan# = 0.577) [a, L33|- The error bars are the 
statistical errors based on the integrated luminosity shown in 
the figure. 



SM. The main challenge would be to reduce the measure- 
ment errors sufficiently to discriminate between models 
and make accurate measurements of the b- and i-quark 
couplings to the Z' ' . The major unknown in the analysis 
is the detection efficiency of the t and &-quarks. To ac- 
count for this we considered two scenarios, an optimistic, 
high efficiency scenario using larger values for e t and e?, 
given in the literature, and a pessimistic, low efficiency 
scenario which used more conservative values. We expect 
that the LHC experiments will attain values somewhere 



in between. Given the promise of this approach, a more 
detailed detector level study to see the effects of detector 
resolution is warranted. 
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